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Abstract 

The white noise approach to the investigation of the dynamics of a quantum particle 
interacting with a dilute and in general non-equilibrium gaseous environment in the 
low density limit is outlined. The low density limit is the kinetic Markovian regime 
when only pair collisions (i.e., collisions of the test particle with one particle of the 
gas at one time moment) contribute to the dynamics. In the white noise approach 
one first proves that the appropriate operators describing the gas converge in the 
sense of appropriate matrix elements to certain operators of quantum white noise. 
Then these white noise operators are used to derive quantum white noise and quantum 
stochastic equations describing the approximate dynamics of the total system consisting 
of the particle and the gas. The derivation is given ab initio, starting from the exact 
microscopic quantum dynamics. The limiting dynamics is described by a quantum 
stochastic equation driven by a quantum Poisson process. This equation then applied 
to the derivation of quantum Langevin equation and linear Boltzmann equation for the 
reduced density matrix of the test particle. The first part of the paper describes the 
approach which was developed by L. Accardi, I.V. Volovich and the author and uses 
the Fock-antiFock (or GNS) representation for the OCR algebra of the gas. The second 
part presents the approach to the derivation of the limiting equations directly in terms 
of the correlation functions, without use of the Fock-antiFock representation. This 
approach simplifies the derivation and allows to express the strength of the quantum 
number process directly in terms of the one-particle 5-matrix. 

*This preprint is a minor modification of the paper published in QP-PQ: Quantum Probability and White 
Noise Analysis 18, Eds. M. Schiirmann & U. Franz, (2005) 428-447. The Appendix with several comments 
is added. The paper was written while the author visited the Centre Vito Volterra of Rome University 
Tor Vergata, Rome, Italy in February-March 2004 on leave from Steklov Mathematical Institute of Russian 
Academy of Sciences. 
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1 Introduction 



The fundamental equations in quantum theory are the Heisenberg and Schrodinger 
equations. However, it is a very difficult problem to solve explicitly these equations for 
realistic physical models and one uses various approximations or limiting procedures 
such as weak coupling and low density limits. These scaling limits describe the long 
time behavior of physical systems in different physical regimes. 

One of the methods to study the long time behavior in quantum theory is the 
stochastic limit method, which was developed by Accardi, Lu and Volovich in the 
book The book is devoted mainly to the weak coupling regime, where one considers 
the long time dynamics of a quantum open system weakly interacting with a reservoir. 
The dynamics of the total system in this limit is described by the solution of a quantum 
stochastic differential equation driven by a quantum Brownian motion; the reduced 
dynamics of the system is described by a quantum Markovian master equation. Also in 
the low density regime the dynamics is described by a quantum stochastic differential 
equation. In this regime one considers the long time quantum dynamics of a test 
particle interacting with a dilute gas (we consider Bose gas in this paper) in the case 
the interaction is not weak but the density n of particles of the gas is small. The 
low density limit is the limit as n ^ 0, t ^ oo, nt = Const {t denotes time) of 
certain quantities such as matrix elements of the evolution operator or the reduced 
density matrix. In this limit the reduced time evolution for the test particle will be 
Markovian, since the characteristic time ts for appreciable action of the gas on the 
test particle (time between collisions) is much larger than the characteristic time tR 
for relaxation of correlations in the gas. 

Accardi and Lu |2l E] and later Rudnicki, Alicki, and Sadowski proved that 
the matrix elements of the evolution operator in the "collective vectors" converge in 
the low density limit to matrix elements of the solution of a quantum stochastic differ- 
ential equation driven by a quantum Poisson process. The quantum Poisson process, 
introduced by Hudson and Parthasarathy [S], should arise naturally in the low density 
limit, as conjectured by Frigerio and Maassen jjj. 

The stochastic golden rule for the low density limit, that is a set of simple rules 
for the derivation of the limiting quantum white noise and stochastic differential equa- 
tions, was developed in Refs [HI El E|) where the case of a discrete spectrum of the 
free test particle's Hamiltonian was considered. This method is based on the white 
noise approach and uses the stochastic limit technique 1^ (in Ref. [TJ| the method is 
generalized to the case of a continuous spectrum). As the main result, the normally 
ordered quantum white noise and stochastic differential equations were derived. Then 
these equations were applied to the derivation of the quantum Langevin and master 
equations describing the evolution of the test particle. 

The idea of the white noise approach is based on the fact that the rescaled free 
evolution of appropriate operators of the gas converges, in the sense of convergence of 
correlation functions, to certain operators of quantum white noise. Then, using these 
limiting operators which are called master fields, one can derive the quantum white 
noise equation for the limiting dynamics and put this equation to the normally ordered 
form which is equivalent to a quantum stochastic differential equation. 

The white noise approach to the derivation of the stochastic equations in Refs. HUHl 
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1^ IHl El uses the Fock-antiFock representation for the canonical commutation relations 
(CCR) algebra of the Bose gas, which is unitary equivalent to Gel'fand-Naimark-Segal 
(GNS) representation. The white noise approach of Ref. jlUj considerably simplifies 
the calculations by avoiding use of the Fock-antiFock representation. A useful tool 
is the energy representation introduced in Refs. |Hlin]) where the case of orthogonal 
formfactors was considered. This consideration was extended in Ref. ^H] to the case of 
arbitrary formfactors and arbitrary, not necessarily equilibrium, quasifree low density 
states of the gas. 

In Ref. |10] to each initial low density state of the Bose gas in the low density 
limit one associates a special "state", which is called a causal state, on the limiting 
master field algebra. The time-ordered (or causal) correlators of the initial Bose field 
converge to the causal correlators of the master fields, which are number operators of 
quantum white noise. This approach allows to express the intensity of the quantum 
Poisson process directly in terms of the one-particle S'-matrix. In this case the algebra 
of the master fields H22|) . the limiting equation (|3U|) . and the quantum Ito table H32|) 
do not depend on the initial state of the Bose gas (see Sect.EJfor details). Instead, the 
information on the initial state of the gas is contained in the limiting state ^pi oi the 
master field [defined by (|25j) - (|2<S|) ]. which is now not the vacuum. Here the operator 
L determines the initial density of particles of the gas; if the gas is at equilibrium 
with inverse temperature (3, then L = e~^^'^ (see the next section). To get the master 
equation one has to take the conditional expectation with respect to the state ^pL. 

The dynamics in the low density limit is given by the solution of the quantum white 
noise equation ()3U() or of the equivalent quantum stochastic equation 

dUt = dNtiS - l)Ut (1) 

where Ut is the evolution operator at time t, S the one-particle S matrix describing 
scattering of the test particle on one particle of the gas, and Nt{S — 1) the quantum 
number process with strength S — 1. In order to describe these objects let us introduce 
two Hilbert spaces TCs and TCi, which are called in this context the system and one- 
particle reservoir Hilbert spaces, and the Fock space r(L^(M_|_; ?^i)) over the Hilbert 
space of square-integrable measurable vector-valued functions from = [0, oo) to 
Til. Then the solution of the equation is a family of operators Ut;t > in 7is 
r(L^(M_l_; T^i)) (adapted process); 5 is a unitary operator in TLs (Xi Hi. 

Let us define the number process. Let X be a self-adjoint operator in a Hilbert space 
/C; for any / G /C let ^{f) be the normalized coherent vector in the Fock space r(/C). 
The number operator N{X) is the generator of the one-parameter unitary group T{e^''^) 
characterized by T{e^''^)"^{f) = "^{e^''^ f); t E M. The number operator is characterized 
by the property {^{f), N{X)^{g)) = {f, Xg){^{f),^{g)) . The definition of N{X) is 
extended by complex linearity to any bounded operator X on /C. Let us consider /C 
of the form L'^{R+;ni) = -L^(M+) ® Tii. For any bounded operators Xq £ B{ns), 
Xi e B{ni) and for any f > define Nt{Xo Xi) := Xq (g) N{x[o,t] ® Xi) and 
extend this definition by linearity to any bounded operator K in 7is ^Tii. The family 
{Nt{K)^t>G of operators in "Hs 'S?'r(L^(M_|_; Hi)) is called quantum number process with 
strength K. 

Equations and H16|) describe the dynamics of the total system and can be applied, 
in particular, to the derivation of the irreversible quantum linear Boltzmann equation 
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for the reduced density matrix of the test particle. This equation can be directly 
obtained from the quantum Langevin equation [U] . 

The reduced dynamics of the test particle in the low density limit with methods, 
based on a quantum Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy, has 
been investigated by Diimcke jl2j . where it is proved that, under some conditions, the 
reduced dynamics is given by a quantum Markovian semigroup. 

In the white noise approach the reduced dynamics can be directly derived from the 
solution of the limiting quantum stochastic differential equation. Namely, the limiting 
evolution operator Ut and the limiting state ipi determine the reduced dynamics by 



where X is any observable of the test particle, (fii') denotes the conditional expec- 
tation, and Tt is the limiting semigroup. This equality shows that Ut is a stochas- 
tic dilation of the limiting Markovian semigroup. Using the quantum Ito table for 
stochastic differential dNt one can derive a quantum Langevin equation for the quan- 
tity U^{X \)Ut- Then, taking partial expectation, one gets an equation for Tt{X) 
and obtains the generator of the semigroup (see the end of Sec. |1J). This is a general 
rule of the white noise approach: one first obtains the Langevin equation and then 
gets the reduced dynamics of the test particle. Let us note that although the quantum 
stochastic equation (fTHj) . which was derived in Refs. 12112], is different from it gives 
the same reduced dynamics. 

The low density limit can be applied to the model of a test particle moving through 
an environment of randomly placed, infinitely heavy scatterers (Lorentz gas) (see the 
review of Spohn^^). In the Boltzmann-Grad limit successive collisions become inde- 
pendent and the averaged over the positions of the scatterers the position and velocity 
distribution of the particle converges to the solution of the linear Boltzmann equation. 
An advantage of the stochastic limit method is that it allows us to derive equations not 
only for averaged over reservoir degrees of freedom dynamics of the test particle but 
for the total system+reservoir. The convergence results and derivation of the linear 
Boltzmann equation for a quantum Lorentz gas in the low density and weak coupling 
limits are presented in Refs. ^^E]- The Coulomb gas at low density is considered in 



The main results presented in this paper are: the causally normally ordered quan- 
tum white noise equations (fT7)|) . and equivalent quantum stochastic equation 
for the limiting evolution operator; the quantum Langevin equation H17|) for the evo- 
lution of any test particle's observable; the linear Boltzmann equation for the reduced 
density matrix (Theorem 

The structure of the paper is the following. In Sec. 2 a test particle interacting with 
a Bose gas is considered. In Sec. 3 the white noise approach developed by L. Accardi, 
I. Volovich, and the author in Refs. 012] is presented. Sec. 4 describes the white noise 
approach developed in Ref. [TU] . 





Ref. dS]. 
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2 Test Particle Interacting with a Dilute Bose 
Gas 



Consider two non-relativistic particles, with masses M and m, which are called the test 
particle and a particle of the gas. Suppose the particles interact by a pair potential 
U{R — r), where R and r denote positions of particles. Then the classical dynamics of 
the particles is determined by the Hamiltonian H^i = P^/2M + p'^/2m + U{R — r), 
where P and p are momentums of the particles. 

The quantum Hamiltonian of such system is obtained by identification of P with 
the momentum operator P, R with the position operator Q and, if instead of one par- 
ticle of mass m there is gas of these particles, by second quantization of the particles of 
the gas. This Hamiltonian has the form H = Hs + Hji + -ffmti where Hs + -ff_R =: Hq 
is the free Hamiltonian with Hs = P'^/2M, Hr = J uj{p)a'^ {p)a{p)dp, uj{p) = p^ /2m] 
the interaction Hamiltonian is H{^t = J U{Q — r)a~^{r)a{r)dr. The free Hamiltonian 
of one particle of the gas Hi is the multiplication operator by the function to. Bo- 
son annihilation and creation operators a{k),a~^(p) satisfy the canonical commutation 
relations 

[aik),a+ip)]=6{k-p) 

The coordinate representation for these operators is introduced as o(r) = J e^^^a{k)dk. 
The interaction Hamiltonian can be written using the Fourier transform of the inter- 
action potential U{p) := J U {r)e^^^ dr as 

Hir,t = j dkdpU{p)e'P^ a+(A; - p)a{k) (3) 

This Hamiltonian acts in the Hilbert space TCs T(TCi), where Tis = Tii = L^(M^). 

Important features of this system are that the interaction Hamiltonian H[^t quadratic 
in creation and annihilation operators and commutes with the number operator, i.e., 
it preserves the number of particles of the gas. 

We will consider, instead of ©, a different interaction Hamiltonian, which nev- 
ertheless keeps its basic properties. More precisely, we consider Hamiltonians of the 
form 

Hir^t A+{go)A{gi) + D+ <g) A+{gi)A{go) 

where D is a bounded operator in Tis; go,gi S Tli are two form-factors, and A{go) = 
f dkgQ{k)a{k) is the smeared annihilation operator. This Hamiltonian is also quadratic 
in creation and annihilation operators and preserves the number of particles of the gas. 
In the present paper we consider the case of a discrete spectrum for the free test 
particle's Hamiltonian, so that 

Hs = Y,enPn (4) 

n 

where e„ is an eigenvalue and P„ is the corresponding projector. This corresponds to 
the situation when the test particle is confined in some spatial region. 

Density of particles of the gas is encoded in the state of the gas, which is chosen to 
be either the Gibbs state at inverse temperature (3, chemical potential fi, and fugacity 
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^ = e^^, or a more general non-equilibrium Gaussian state, i.e., a gauge invariant mean 
zero Gaussian state with the two point correlation function 

V.L,€(^+(/M(5))=e(5,Y^/) (5) 

Here ^ > is a small positive number and L is a bounded positive operator in TCi 
commuting with the one-particle free evolution St = e^^^^ [the multiplication operator 
by a function L{k)]. In the case L = e~^^^, so that L{k) = e~^'^(^), the state (pi^^ is 
just the Gibbs state with the two point correlation function 

<PL,da+{k)a{k')) = n{k)6{k - k') 

where n{k) is the density of particles of the gas with momentum k: 

n{k) ^ 



-I3k^/2m 



1 _ ^g-/3fe2/2m 

Notice that in the limit ^ ^ the density goes to zero. Therefore the limit ^ ^ is 
equivalent to the limit n(fc) 0. 

The dynamics of the total system is determined by the evolution operator which in 
interaction representation has the form U{t) := Q^^^hccQ-^tHtot _ 'pj^g evolution operator 
satisfies the differential equation 

^ = -ii?int(t)f/(t), 

where i?int(i) = e'*''^f"'°<=//inte~'*''^*''°'= is the free evolution of the interaction Hamiltonian. 
The iterated series for the evolution operator is 

^ t tn~l 



(6) 



U{t) = 1 + Y,{-iT I dh... I dtnHi^tih) . . . ^int(t„) 


With the notation D{t) := q^'^^s D^-itHs ^j^g evolved interaction becomes 

Hi^tit) := D{t) ® A+{Stgo)A{Stgi) + D+{t) ® A+ {Stgi)A{Stgo). 

Using the spectral decomposition @ and introducing the set of all Bohr frequencies 
B, that is the spectrum of the free test particle's Liouvillean i[Ha,, •], one can write the 
free evolution of D as 

The reduced dynamics of any test particle's observable X in the low density limit 
is defined as the limit 

Tt{X) := Ihn 99,^,5 (C/+(t/e)(X ® l)U{t/i)) 

where ^PL,i{-) denotes partial expectation. The reduced density matrix p{t) is defined 
through the duality Tr (/>(0)rt(X)) = Tr(/>(t)X). As it was mentioned in the Intro- 
duction, in the white noise approach the generator of the limiting semigroup can be 
easily derived from the quantum white noise equation. 
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3 The White Noise Approach 

In Refs. 121 El m the dynamics of the total system is constructed in the Fock-antiFock 
representation for the CCR algebra, which is unitary equivalent to the GNS represen- 
tation. It is defined as follows. 

Denote by W^. the conjugate of Tii, i.e. TC[ is identified to Tli as a set and the 
identity operator l : TLi — > TLi is antilinear: V/ G Wi, c G C 

.(c/) = c*.(/), {i{f),i{g)). = {g,f) 

Then, Tii is a Hilbert space and, if the vectors of 7ii are thought as ket-vectors |.^), 
then the vectors of Tii can be thought as bra- vectors (^|. The corresponding Fock 
space TiTLi) is called the anti-Fock space. 

In this section we assume that for any t G M: {go,Stgi) = 0. In this case it was 
shown in Ref. 2 that the dynamics of the total system is given by the family of unitary 
operators uj:^^ in TCs r('Hi) (8) r('H^) which satisfy the Schrodinger equation: 

dtup = -iH^mP, = 1- 

Here the part of the modified Hamiltonian, which gives a nontrivial contribution in the 
low density limit, has the form 



+ ^i[A{Stgi) ® AiSte-^^'^'^go) + ^+(5tgo) ® A-^iSte-^^^^'^gi)] } + h.c. 
After the time rescaling t ^ t/^ the evolution operator satisfies the equation 



where we introduce for each n, m = 0, 1 and to ^ B the rescaled fields: 

Nn,m,i{^.t) := -^e-''^I^A+{St,^gn)A{St/^g^)®l (7) 
Bn,m,i{u:,t) := l=e''^I^A{St/^gn)®A{Sti^e-''''^/^gm) (8) 

and is the adjoint of Bn,m,ii^,t)- 

3.1 Master Field 

It is convenient to use the energy representation for the investigation of the limit as 
^ ^ of the rescaled fields (O, It is defined in terms of the projections 

Pe ■■= 5{Hi - E) 
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which satisfy the properties 

PePe' = 5{E - E')Pe, Pe = Pe, St = j dEPEs''^ 
Define the energy representation for the fields Q, © as 

^it{Ei-E2-w)/£, 

Nn,m,^{Ei,E2,UJ,t) := A+{PE,gn)A{PE2gm)®l (9) 

? 

it{E2-E■^+uJ)/^, 

Bn,m,dEi,E2,L0,t) := -j= g„) ® ^(P^^ e"^^'/' <7m) (10) 

and let B^^^{Ei,E2,u},t) be the adjoint of Bn^m,^{Ei, E2,io,t). The operators 0,® 
can be expressed in terms of ©,(^01) as 

Nn,m,ii^^,t) = j dEidE2Nn^m,i{El, E2,UJ,t) 

Bn,m,d'^,t) = j dEidE2B {Ei,E2,uj,t) 

Let be the set of all linear combinations of the Bohr frequencies with integer 
coefficients: Q = {uj\uj = '^^.rikOOk withn^ G Z, G B}. Extend the definition of the 
fields ((ni),(ITni) to arbitrary a; G 17. The limit as ^ ^ of the fields (l9ll. (frnj) was found 
in Ref. <^ and is given by the following theorem. 

Theorem 1 The limits of the rescaled fields 

Xn,miEi,E2,UJ,t) ■.= \imXn,m,dEi,E2,uj,t), X = B , B+ , N 

exist in the sense of convergence of correlators and satisfy the commutation relations 

[Bn,miEi, E2,uj,t), B:^,^^,{Es, E4„uj' ,t')] = 27r(5^_^/ J„^„/(5„_^/(5(t' - t) 
x6{Ei - E^)5{E2 - E4)d{Ei - E2 - uj){gn, PE,gn){gm, PE^e-^^'g^) (11) 

[B^^m{El,E2,iO, t), iV„',^' (^3, ^4, t')] = 2^5^y6{t' - t) 
x5{Ei-E:i)5{Ei-E2-Uj){gn,PE,gn)B^,^m{EA,E2,LO-Uj\t) (12) 
[Nn,UEl,E2,UJ, t),Nn',m'{E3, Ei,Uj\ t')] = 2^6{t' - t) 
X{'^m,n"5(-E2 - E2,)5{E2, - Ei + Lj) {gm, PEigm) N^^rn' {El, E^, U + Uj' , t) 
-6n,m'S{Ei - E4)6{Es - El - Uj'){grr,PE,gn)Nn',m{E3,E2,UJ + Uj',t)} (13) 

The causal commutation relations of the master field are obtained replacing in Ul\) - 
lilS\) the factor 5{t' — t) by 6^{t' — t), where the causal 5 -function 6^{t' — t) is defined in 
Ref ty/, section (8.4); 27r6{Ei- E2-UJ) by {i{Ei- E2-UJ -iO)y^ and 27r6{E3- Ei±lo) 
by {i{E3-Ei±uj-iO))-K 
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3.2 Fock Representation of the Master Field 

The next step is to reahze the algebra of the master field by operators acting in a 
Hilbert space. We realize it in the Fock space, which is constructed as follows. 

Let K he a. vector space of finite rank operators acting on the one-particle Hilbert 
space Til with the property that for any lo £ ^l; X,Y £ K 



{X,Y)^ : = j dtTr (^e~^"'X*StYS;^ e"*"^* = j dtTr (^e-'^^'X*Yt 

= 27r y dETr {e'^"^ X* PeY Pe-^^ < oo 

where Yt = StYS^ the free evolution of Y. It was shown in Ref. that the space 
K is non empty and (•, •)a; defines a prescalar product on K. Let be the Hilbert 
space with inner product obtained as completion of the quotient of K by the 

zero (•,-)a;-norm elements. Denote /C := K^. 

Consider the Fock space 



r(L^(M+) ® /C) = r((^ L^(M+, K^)) = T{L'{R+, K^)) (14) 

where the last infinite tensor product is referred to the vacuum vectors. Let bf^^{X), bt^u){X) 
be the white noise creation and annihilation operators in this Fock space. These oper- 
ators satisfy the commutation relations 

[bt,UX),b+^,{Y)] = 6{t' - t)6^,^,{X,Y)^. 

Each white noise operator bt^uj{') acts as usual annihilation operator in r{L'^{^^, K^)) 
and as identity operator in other subspaces. 

The representation of the algebra ((TT]) ~ H13|) can be constructed in the Fock space (|14() 
by the identification 

Bn,m{El,E2,iO,t) = bt,uj{\PEign){PE2 9m\) 

The number operator is defined as 

Xn,m{El^E2,tO,t) = 

E E li,{Er-u:r)bl^^{\gn){PE, —uj-i_9e\) 
e=0,l wieSl 

with He{E) '■= {9e, Pe^'^^^ Qe)'^^ ■ One easily checks that these operators satisfy the 
commutation relations (|11 |) - (|13|) . Denote 

iV„,„(a;,t) := j dEidE2Nn,m{Ei,E2,uj,t) 

Bn,miE,UJ,t) := j dE'Bn,miE',E,UJ,t) 

The limiting white noise Hamiltonian acts in TLs r(L^(M+) (g) /C) as 

H{t) = ^D^®\^No^i{uj,t)+ I dE^Bi^oiE,-LO,t) + B^^^{E,u;,t)j^ + h.c. 
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3.3 Quantum White Noise, Langevin and Boltzmann Equa- 
tions 

The white noise Schrodinger equation for the evolution operator in the low density 
limit is 

dtUt = -iH{t)Ut (15) 

Following the general theory of white noise equations, in order to give a precise meaning 
to this equation we will put it in the causally normally ordered form, in which all 
annihilators are put on the right hand side of the evolution operator and all creators 
are on the left hand side. This procedure gives a normally ordered quantum white 
noise equation, which is equivalent to a quantum stochastic differential equation. 

For any n,m G {0, 1}, lo £ Cl, E £ M+ let R^'^,{E) be operators in TCs which are 
explicitly defined in Sec. 7, Ref. iSj. The following theorem was proved in Ref. [2]. 

Theorem 2 The normally ordered form of equation ij,5)) is 

dtUt= E IdE[ZK'!i:'iE) E f^e{E)B+,{E,u;,t)UtBmAE,uj',t) 

n,m=0,l ll1,uj' £=0,1 

+ E (K';^{E)B+^{E, u;, t)Ut + Rl^:::{E)UtBn,m{E, ^, t)) 

+R'^;^{E) {gn , Pi^ e-/5^i g^)Ut] (16) 

The normally ordered equation H16() can be applied to derivation of the quantum 
Langevin equation for test particle's observables. Let X be any observable of the test 
particle. The Langevin equation is the equation satisfied by the stochastic flow jt 
defined by jt{X) = Xt := XUf 

Theorem 3 The quantity Xt satisfies the quantum Langevin equation: 

Xt= E !dE\ T.l^e{E)B+^{E,uji,t)UtQZ':%2{X)UtB„,,,{E,uj2,t) 

n,m=0,l ^Lci,LU2 £ 

+ E (B+n,{E, uj, t)C/;G::'J(X)^, + Ut*e^^i^{X)UtB^,^{E, u, t)) 

+(5n,Pi^e-^^^<7m)t/t*e^;o"(X)^t] (17) 
where the structure maps are @Z'^c^^{X) := X R^'^^2(E) + Ki^^i]Jl{E)X + 2'YRe'y£{E + 

e,Lu 

w)i?i',t^"(-E)Xi?5',™2 (E) ■ 

Equation p7j) can be written in terms of the stochastic differentials: 

d3t{X)=jtoY.jdE\ E Ql'^u.2{X)dNt{Z::,C^,{E)) 



n,m 



+E e::;™(x)di?+((i5„)(Pi.5^i).)+G^:rwrf^t((i<7n)(PE<7™i). 

+jt o C{X)dt 
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where Z"j^2(-E') is a certain operator in /C, which is expHcitly defined in Ref. [21; C 
is a quantum Markovian generator, which has the form of a generator of a quantum 
dynamical semigroup [T7|: 

LiX) = - ^{^1/(1), X} + i[Feff, X] 

Here = 2^ ^ / '^^(5., ^i?e-'«^i^/,)(5.'P£+.<7e')<o''(^)^<o (^) is a com- 

pletely positive map and the effective Hamiltonian i^gff := Xle / dE{gi;, PEe~^^^ Qs) x 
- i?o;o(£;))/2i is selfadjoint. 
The following theorem is a direct consequence of the quantum Langevin equa- 
tion (dZI of Ref. % 

Theorem 4 The reduced density matrix satisfies the quantum linear Boltzmann equa- 
tion 

^ = Up{t)) (18) 

where the generator is the dual to the generator C. 

The generator is the sum of its dissipative and Hamiltonian parts, C^:{p) = 
'Cdiss(p) ~ i[HcS,p]- The explicit form can be obtained by direct calculations using the 
expression above for the generator C and Eq. (8.6) of Ref 9 for the T operator. Finally 
it has the following form. Let T be the one-particle T-operator for the scattering of 
the test particle and one particle of the gas and Tn^n'{k,k') = {n, k\T\n' , k') be its 
generic matrix element, where |n) is an eigenvector of Hs with eigenvalue and k the 
momentum of one particle of the gas. Denote 

T^{k,k'):= ^ Tm,n{k,k')\m){n\ 

m,n: em—en=oj 

Density of particles of the gas is determined by the function L{k) in For the Gibbs 
state L(k) = e~^^^^\ Other forms of the density correspond to non-equilibrium states 
of the gas and can be controlled, for example, by filtering. In these notations the 
dissipative part of the generator is 

£diss(/o) = 2^ E / dkdk'5{u{k') - uj{k)+uj)L{k) 
ujeB 

TUk',k)pT+{k',k) -^(T+ik',k)TUk',k)p + pT+ik',k)TUk',k)^] (19) 

In the case the gas is in equilibrium this generator coincides with the generator of the 
quantum linear Boltzmann equation obtained in Ref. jl2j . 

4 White Noise Approach without Fock-antiFock 
Representation 

The approach to derivation of the quantum white noise equations directly in terms 
of the correlation functions, without use of the Fock-antiFock representation, was de- 
veloped in Ref. [TU!. In this approach one introduces the notion of causal state and 
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causal time-energy quantum white noise and proves the convergence of chronological 
correlation functions of operators 

Nf,g,^{t) = ^A+{St/^f)A{St/^g) (20) 

acting in T{TCi) to correlation functions of the time-energy quantum white noise (The- 
orem [SJ. This time-energy quantum white noise is a family of creation and annihilation 
operators, with commutator proportional to (5-function of time and energy [see H22|l ]. 
These operators act in a Fock space which, in difference with (fTl|) . does not depend on 
the initial state of the gas 

Suppose for simplicity that D{t) = D. Then the evolution operator U{t/^) after 
the time rescaling t ^ t/^ satisfies the equation 

= _,(Z) iV,„,,„5(t) + D+^ N,,,g,,^it))Um) (21) 

Theorem 1^1 and the causal commutation relations (|2c{|) are used to show that the limit 
as ^ — > of the rescaled evolution operator satisfies the causally normally ordered 
equation (|30|) . That equation is equivalent to the quantum stochastic equation (|3T|) 
which can be written in Hilbert module notations as (|34() and then in terms of the 
one-particle S-matrix as (|36|) . 

As it was stated in the Introduction, in this approach the algebra of the time-energy 
quantum white noise (|22|) . the quantum Ito table p2|) and the quantum stochastic 
equation for the limiting evolution (|36jl do not depend on the initial state of the gas. 
This is different from the approach of Sect. |31 where the commutation relations 
for the master field, the Hilbert space representation (subsection 13. 2 j) and hence the 
limiting equation (|16() depend on the initial state (through the factor e~^^^). Instead, 
the dependence on the initial state of the gas now is contained in the limiting state 93^ 
(in the equilibrium L = e~^^^) (in the approach of Sect. |31 the state of the master field 
is the vacuum state). Considering the limiting state as the conditional expectation 
[with the property H33() ]. one can derive the quantum master equation for the reduced 
dynamics of the test particle which coincides, when restricted to the same model, with 
the analogous equation following from H17|) . 



4.1 Causal Time-Energy Quantum White Noise 

Define the Hilbert space Xt-^^^^Hi as the completion of the quotient of the set 

|F:M+^?^l s.t. :=27r j dE{F{E), PeF{E)) < 00^ 

with respect to the zero-norm elements. The inner product in ^^7^1,//! is {F, G) = 
27r / (1E{F{E),PeG{E)). Let Bj{E,t), Bg{E',t') be the creation and annihilation 
operators acting in the symmetric Fock space r(L^(M+, A'^^^^^/f^)) over the Hilbert space 
L^(M+, AV^i,Hi) of square integrable functions / : M+ — > Xj^-^^^h-^. These operators 
(operator-valued distributions) satisfy the canonical commutation relations 

[Bg{E,t), B+{E',t')] = 27T6{t' -t)6{E' - E){g,PEf) (22) 
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and causal commutation relations 



[B,{E,t), B+{E' ,t')] = 5+{t' -t)5{E' - E)^gj{E) (23) 

where 5+(t' — t) is the causal (5-function and jgj{E) = f^^dt{g, Stf)e~''^^ . The 
meaning of two different commutators (|22|) and (|23|) for the same operators is explained 
in Ref. [l], Sect. 7. These operators are called time-energy quantum white noise due to 
the presence of 6{t' - t)6{E' - E) in (|^ . 
Define the white noise number operators as 

Nf,g{t) = J dEB+{E,t)Bg{E,t) (24) 

For any positive bounded operator L in TCi define the causal gauge-invariant mean- 
zero Gaussian state by the properties (|25|) - l|28|) : 

for n = 2k ^LiB^...B^)=Y,VLiB:^B;^)...MB::'B':^) (25) 

where the sum is taken over all permutations of the set (1, . . . , 2k) such that ia < ja, 
a = 1, . . . ,k and ii < i2 < • • • < ik] B^ := {Em, tm) for m = 1, . . . , n are time- 
energy quantum white noise operators with causal commutation relations (|23p . and c^, 
means either creation or annihilation operator; 

forn = 2k + l v?l(5^i . . . fi^") = (26) 



<fL{Bf{E,t)Bg{E',t')) = VL{Bj{E,t)B+{E',t')) = (27) 

^L{B^{E,t)Bg{E,t')) = Xm(t'){9,PELf) (28) 

The "state" ipi does not satisfy the positivity condition. This is a well-known situation 
for the weak coupling limit (see Ref. (Ij ) and is due to the fact that we work with time- 
ordered, or causal correlators. Therefore it is natural to call such "states" causal states. 

Def. 1 Causal time-energy quantum white noise is a pair {B^{E, t), (pi), where B^{E, t) 
satisfy the causal commutation relations \2'J\) and if l is a causal gauge-invariant mean- 
zero Gaussian state. 

Theorem 5 For any n G N in the sense of distributions over simplex ti > t2 > ■ ■ ■ > 
^ one has the limit 

hm^L,dNf,^9idtl) ■ ■ ■ Nf„,g„dtn)) = ^l{N f,,g,{ti) . . . iV/„,g„(t„)) 

This theorem was proved in Ref. jlUj . 

Remark 1 This convergence is called convergence in the sense of time-ordered cor- 
relators. The fact that we use the distributions over simplex is motivated by iterated 
series ^ for the evolution operator. 
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Remark 2 The proof is based on the fact that for any n £ N only one connected 
diagram survives in the limit. This can be interpreted as emergence of a new statistics 
( different from Bose ) in the low density limit^ . For a discussion of new statistic arising 
in the weak coupling limit see Ref. 

The following theorem is important for investigation of the limiting white noise 
equation for the evolution operator. 

Theorem 6 The limit state ipi has the following factorization property: Vn G N 

{Bf{E,t)Nf,^g,{t,) . . . Nf^^g^{tn)Bg{E,t)) 

= ipL{Bl{E,t)Bg{E,t))^L{Nf,^gAti) . ■ ■ Nf^,gM) 

where the equality is understood in the sense of distributions over simplex t > ti > 
t2 > • • • > in > 0. 

Theorem IS] allows us to calculate the partial expectation of the evolution operator 
and Heisenberg evolution of any system observable in the low density limit. In fact, 
partial expectation of the n-th term of the iterated series © (or equivalent series for 
Heisenberg evolution of a system observable) after time rescaling t ^ t/^ includes the 
quantity 

t tn-l 

dtl... J dtn^L,d^f,,,,,iitl) . . . Nf„^g„^^itn)) 


(where fa, ga are equal to go or gi). The limit as ^ of this quantity can be 
calculated using Theorem El For example, the contribution of the connected diagram 
is equal to 

j dti j dt25+{t2-ti) j dt36+it3 - t2) . . . j dtn5+{tn-tn-l) 


X j dE{gn, PELfi)jg,j,{E) . . . jg„_„u{E) 
= t I dE{gn,PELfihg,j,{E)...la„,ijSE) 

Similarly one can calculate the contribution of nonconnected diagrams (they give terms 
proportional to higher powers of t). Summation over all orders of the iterated series 
gives the reduced dynamics of the test particle. An advantage of the white noise 
approach is that it allows to get the limiting dynamics in a nonperturbative way, 
without direct summation of the iterated series. This procedure includes derivation of 
the causally normally ordered white noise equation for the limiting evolution operator. 
After that the reduced dynamics of the test particle can be easily found. 



^In Ref. |18| the statistics which appears in the low density limit was found and a connection with 
Voiculescu free independence theory was established. 
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4.2 The White Noise and Quantum Stochastic Equations 

The limiting evolution operator satisfies the white noise Schrodinger equation 

^ = -i{D iVgo,9i W + Ng,,g,{t))Ut, Uo = 1 (29) 

The next step is to put this equation to the causally normally ordered form, i.e., to 
put all annihilation operators, appearing in Nf^g{t), on the right side of the evolution 
operator and all creation operators on the left side. 

Assume that for each E £M the following inverse operators exist 

To{E) := (l + 790,51 (^)^'' - lgu9oiE)D + {lg,,g,lg,,g, - lgu9o79o,9i)iE)DD+y' 

-1 



Tl{E) := ^1 + 7go,gAE)D+ - lg„g,{E)D + {lg,,g,lg„g, " ^g,,g,^g,,g,){E)D+ D 

Denote 

R^fl{E) := jg,,g,{E)DTi{E)D+, Ro,i{E) := -DTi{E){l + jg,,g,{E)D- 
Ri,i{E) := ^g„g,{E)D+n{E)D, R,,o{E) := D+To{E){l - ^g,,g,{E)D) 

Theorem 7 The causally normally ordered form of equation \29\) is 



dt 



I dER,n,n{E)B^jE,t)UtBg^{E,t) 

n,m=0,l 



(30) 



This theorem was proved in Ref. jlUj . 



Remark 3 An immediate consequence of Theorem is the following factorization 
property of the limiting state ipL-' 

ipL{B+{E,t)UtBg{E,t)) = ipL{B+{E,t)Bg{E,t))ipL{Ut) 

This property of the state ip^ similar to the factorization property of the state deter- 
mined by a coherent vector ^, ||^|| = 1.' 

{^,B+{E,t)UtBg{E,t)^) = {^f,B+{E,t)Bg{E,t)^){^,Ut^) 

which is usually used to define quantum stochastic differential equations. 

Normally ordered white noise equation (|3Ujl equivalent through identification 

B+iE,t)UtBniE,t)dt = 27rdNti\PEgm){PEgn\)Ut 
to the quantum stochastic differential equation 

dUt = -2TT f dERmAE)dNt{\PEgm){PE9n\)Ut (31) 

n,m=0,l 

where Nt is the quantum number process in T{L^{R+)iE)ni). The stochastic differential 
dNt satisfies the quantum Ito table 

dNtiX)dNt{Y) = dNtiXY) (32) 
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where X,Y are operators in TCi. The hmiting state ipi has the property 

ipL{27rdNt{\PEf){PEg\)) = {9,PELf)dt (33) 
Equation (|31() can be written in Hilbert module notations as 

dUt = dNt(^-27r j dERm,n{E)(^\PE9m){PEgn\)Ut (34) 

n,m=0,l 

The one-particle 5-matrix for scattering of the test particle on one particle of the gas 
has the form 

5 = l-27r j '^ERrn,n{E)®\PEgm){PE9n\ (35) 

n,7Tt=0,l 

This is a unitary operator: S = SS^ = 1. An immediate conclusion from ()34() and 
(|35)) is the following theorem which was proved in Ref. JHl and is one of the main 
results of the paper. 

Theorem 8 The evolution operator in the low density limit satisfies the quantum 
stochastic equation driven by the quantum number process with strength S — 1: 

dUt = dNtiS - l)Ut (36) 

This equation can be applied to derivation of the quantum master equation for the 
reduced dynamics. Let X £ B{Ti.s) be an observable of the test particle. Its time 
evolution Xt = XUt satisfies the equation 

dXt = dU^XUt + U^XdUt + dU^XdUt 

(we identify X and X Now, since 

dC/+ = U^dNt{S+ - 1) 

(we use the Hilbert module notations, hence dNt{S^ — 1) does not commute with U^) 
and using the quantum Ito table (|32|) one gets 

dXt = U+dNt{S+ -l)XUt + U+XdNt{S -l)Ut 

+U+dNtiS+ - l)XdNtiS - l)Ut = U+dNt{e{X))Ut 

where the map : S(?^s) ^ 5(Ws ® Wi) has the form Q{X) = {S+ - l)X{S - 1) + 
(5"*" — 1)X + X{S — 1) = S^XS — X. Simple computations, together with the explicit 
form H35|l for the S-matrix and the property PePe' = ^{E — E')Pe of the projector, 
give the expression 

e(X) = 27r J2 I dEQ^'jf'iX) ® \PE9^m){PEgn\ 

n,m=0,l 

where G^'™(X) = 27r E R^, ^JE)XR^> ^n{E){9n'PEgm') - Rlm{E)X - XRm,n{E). 

n' ,m' 

Therefore 

dXt = 27: Yl j dEdNt{\PEgm){PE9n\)UtQT'{X)Ut 
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The reduced dynamics Xt := (piiXt) G B{7is) is obtained by taking conditional 
expectation of both sides of this equation in the state ipi and using 



When restricted to the case of orthogonal form- factors go,gi, this master equation 
coincides with the one of Ref. 0. The linear Boltzmann equation for the reduced 
density matrix is the dual to this equation and has the form of Eq. H18() with dissipative 
generator of the form Eq. (|19() . 



This paper presents two versions of the white noise approach to the investigation of the 
dynamics of a quantum system interacting with a gaseous environment. The approach 
allows to describe the dynamics of the total system consisting of the particle and the 
environment in the low density limit. The low density limit corresponds to the kinetic 
regime when only pair collisions (i.e., collisions of the test particle at one time moment 
with one particle of the gas) contribute to the dynamics such that the probabilities 
of multi particle collisions are negligible. Although the main result of the approach is 
the quantum stochastic differential equation for the total dynamics, quantum master 
equation for the reduced system dynamics is derived as its simple consequence [Eq. (fT8|) 
with Lindblad-GKS generator given by Eq. (|19() ]. The derivation is performed ab initio 
for general, including non-equilibrium, gases. 

Section 2 of the present paper considers cases of discrete and continuous spectrum 
of the test particle. The limiting equations are explicitly derived here for the case of 
discrete spectrum. This corresponds to confining the test particle in a spatially finite 
region, which however can be arbitrarily large. The master equation (|18j) describing 
the dynamics of a particle interacting with a gas recently was applied to the problem of 
non- unitary quantum control U^] j where it was used to analyze capabilities of control- 
ling quantum systems by optimizing with learning control algorithms the state of the 
surrounding gas (i.e., by optimizing its distribution function). Master equations for 
a quantum Brownian particle in a free space (i.e., with translation invariant dynam- 
ics) also attract attention of the researchers. In this context I mention recent works 
by Bassano Vacchini 120], Klaus Hornberger [2^, and Stephen Adler [22j, where also 
other relevant references can be found. The two characteristic features of the white 
noise approach are that it allows to derive solvable equations for the dynamics of the 
total system (the test particle and the environment) and that these equations are not 
phenomenological, they are derived from the exact microscopic dynamics. 
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